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CRITICAL  COMBINATIONS  OF  SHEAR  AND  TRANSVERSE  DIRECT  STRESS  FOR  AN  INFINITELY 
LONG  FLAT  PLATE  WITH  EDGES  ‘ELASTICALLY  RESTRAINED  AGAINST  ROTATION 

By  S.  B.  BA.TDORF  and  John  C.  HooBotT 


SUMMARY 

An  exact  solution  and  a  closely  concurring  approximate 
energy  solution  are  given  for  the  buckling  of  an  infinitely  long 
flat  plate  under  combined  shear  and  transverse  direct  stress 
with  edges  elastically  restrained  against  rotation.  It  was  found 
that  an  appreciable  fraction  of  the  critical  stress  in  pure  shear 
may  he  applied  to  the  plate  without  any  reduction  in  the  trans¬ 
verse  compressive  stress  necessary  to  produce  buckling.  An 
interaction  formula  in  general  use  was  shown  to  be  decidedly 
conservative  for  the  range  in  which  it  is  supposed  to  apply. 

INTRODUCTION 

In  the  design  of  stressed-skin  structures,  consideration 
must  sometimes  be  given  to  tbo  critical  stresses  for  a  sheet 
under  a  combination  of  shear  and  direct  stress.  The  upper 
surface  of  a  wing  in  normal  flight,  for  example,  is  subjected 
to  combined  shear  and  compressive  stress,  and  the  lower 
surface  is  subjected  to  combined  shear  and  tensile  stress. 
The  upper  surface  may  then  buckle  at  a  lower  compressive 
or  shear  stress  than  if  either  stress  were  acting  alone.  The 
critical  shear  stress  for  the  lower  surface  may  be  increased 
by  the  presence  of  the  tensile  stress. 

If  the  wir^  has  closely  spaced  chordwise  stiflEeners,  the 
skin  between  two  adjacent  stififeners  may  be  regarded  as  a 
long  sheet  slightly  curved  in  the  longitudinal  or  chordwise 
direction,  straight  in  the  transverse  or  spanwise  direction, 
and  loaded  in  shear  and  transverse  direct  stress.  A  con¬ 
servative  preliminary  estimate  of  the  critical  stresses  may  be 
obtained  if  the  sheet  is  considered  to  be  flat  and  infinitely 
long.  In  the  present  paper,  the  critical  stresses  are  computed 
for  an  infinitely  long  flat  plate  loaded  as  indicated  in  figure 
1  (a).  The  corresponding  idealization  of  the  case  of  a  wing 
with  spanwise  stiffeners,  sho^vn  in  figure  1  (b),  was  treated 
in  reference  1. 

CONVENTIONAL  INTERACTION  FORMULAS 

For  buckling  of  structures  under  combined  loading  con¬ 
ditions,  no  general  theory  has  been  developed  that  is  appli¬ 
cable  to  aU  cases.  Stress  ratios  (reference  2),  however, 
provide  a  convenient  method  of  representing  such  conditions. 
For  example,  the  ratio  of  the  shear  stress  actually  present 
in  a  structure  to  the  critical  shear  stress  of  the  structure 
when  no  other  stresses  are  present  may  be  called  the  shear- 
stress  ratio.  Stress  ratios  may  aimilarly  be  defined  for  each 
type  of  stress  occurring  in  the  structure. 
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Ca)  Type  of  loading  problem  solved  In  present  paper. 

(b)  Type  of  loading  problem  solved  In  reference  1. 

FIOCEE 1.— Buckling  of  an  InOnltely  long  plate  under  combined  loads. 

It  is  generally  assumed  that  equations  of  the  type 

...  =1  (1) 

may  be  used  to  express  tbe  buckling  conditions  in  the  case 
of  combined  loading  (reference  3,  pp.  1-18) .  In  equation  (1) , 
R\,  i?2,  and  B*  are  stress  ratios  and  p,  g,  and  r  are  exponents 
chosen  to  fit  the  known  results.  (All  s3fmbols  are  defined  in 
appendix  A.)  Such  a  formula  gives  the  correct  results  when 
only  one  type  of  loading  is  present  and  has  the  further 
advantage  of  being  nondimensional.  Equation  (1)  implies, 
moreover,  that  the  presence  of  any  positive  fraction  of  the 
critical  stress  of  one  type  reduces  the  amount  of  another 
type  of  stress  required  to  produce  buckles;  this  implication 
appears  reasonable  and  has  been  proved  tiue  in  some  cases 
(references  1,  2,  and  4). 

In  reference  2  the  following  interaction  formula  is  given 
for  an  mfinitely  long  plate  with  clamped  edges  loaded  in 
shear  and  longitudinal  compression: 

R,^-^+R,=  l  (2) 

where  R,  is  shear-stress  ratio  and  Re  is  longitudinal  direct- 
stress  ratio.  The  same  formula  is  recommended  in  reference 
3  for  general  use  for  the  buckling  of  any  flat  rectangular 
plate,  regardless  of  the  direction  of  compression  and  the 
degree  of  edge  restraint. 

Later  theoretical  work  (reference  1)  shows  that,  to  a  high 
degree  of  accuracy,  for  an  infinitely  long  plate  with  any 
degree  of  edge  restraint  loaded  in  shear  and  longitudinal 
compression 

R,^+Re=l  (3) 
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Fioure  2.— Crltioal  combinations  ot  sbear.stress  coefficient  k,  and  transveise  direct-stress  coefficient  for  an  infinitely  long  plate  with  edges  elastically  restrained  against  rotation. 
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The  same  foiinula  was  found  in  reference  5  to  be  applicable 
to  simply  supported  rectangular  plates  of  aspect  ratios 
0.5,  1,  and  2;  the  conclusion  was  drawn  that  interaction 
curves  in  sti’css-ratio  form  are  practically  independent  of 
the  dimensions  of  the  plate. 

The  present  analysis,  however,  indicates  that  the  buckling 
of  an  inlinitoly  long  plate  loaded  in  shear  and  transverse 
compression  (fig.  1  (a))  is  not  adequately  represented  either 
by  equation  (2)  or  (3)  or  by  any  formula  of  the  type  of 
equation  (1).  Two  independent  theoretical  solutions  to 
this  buckling  problem  are  given  in  appendixes  B  and  C. 
Appendix  B  contains  the  exact  solution  of  the  differential 
equation  of  equilibrium,  and  appendix  C  contains  an  ePeigy 
solution  leading  directly  to  an  interaction  foimula.  This 
energy  solution,  which  gives  approximate  values  only,  was 
made  to  obtain  an  initial  quick  survey  of  the  problem  and 
to  provide  a  check  on  the  results  of  the  exact  solution. 
Approximate  interaction  formulas  in  substantial  agreement 
with  these  results  were  given  for  the  cases  of  simply  sup¬ 
ported  and  clamped  edges  in  reference  6. 

RESULTS  AND  DISCUSSION 

In  figure  2,  curves  are  given  that  indicate  the  critical 
combinations  of  shear  and  transverse  direct  stress  for  an 
infinitely  long  plate  with  edges  elastically  restrained  against 
rotation.  These  curves  are  computed  from  the  exact  solu¬ 
tion  presented  in  appendix  B.  The  degree  of  edge  restraint 
is  denoted  by  e,  which  is  defined  in  appendix  B  in  such  a 
way  that  zero  edge  restraint  corresponds  to  simply  supported 
edges  and  infinite  edge  restraint  indicates  clamped  edges. 
A  similar  set  of  curves  is  given  in  terms  of  stress  ratios  in 
figure  3.  The  numerical  values  used  to  plot- figures  2  and  3, 
together  with  the  values  found  by  the  energy  solution,  are 
given  in  table  I. 

The  most  striking  feature  of  these  results  is  that  an  appre¬ 
ciable  fraction  of  the  critical  stress  in  pure  shear  can  evi¬ 


dently  be  applied  to  the  plate  without  any  reduction  in  tiie 
compressive  stress  necessary  to  produce  buckling.  (See 
fig.  3.)  This  fraction  varies  from  about  one-third  to  more 
than  one-half,  depending  on  the  degree  of  restraint.  At 
shear  stresses  higher  than  those  corresponding  to  this  frac¬ 
tion,  the  compressive  stress  required  to  produce  buckling  is 
reduced  by  the  presence  of  shear.  The  result  that  the  com¬ 
pressive  buckling  stress  is  entirely  unaffected  by  the  presence 
of  a  considerable  amount  of  shear  is  probably  peculiar  to 
infinitely  long  plates.  It  is  to  be  expected,  however,  that 
this  result  will  be  closely  approached  in  the  case  of  long 
finite  plates. 

In  figure  4  a  comparison  is  made  between  the  exact  solu¬ 
tions  "and  the  interaction  formulas  of  equations  (2)  and  (3). 
Equation  (2),  which  is  the  interaction  formula  in  general 
use,  is  seen  to  be  decidedly  consei’vative. 

CONCLUSIONS 

The  exact  solution  of  the  differential  equation  for  the 
buckling  of  an  infinitely  long  flat  plate  under  combined  shear 
and  transverse  direct  stress  with  edges  elastically  restrained 
against  rotation  indicates  the  following; 

1.  An  infinitely  long  flat~plate  may  be  loaded  with  an 
appreciable  fraction  of  its  critical  stress  in  pure  shear  without 
causing  any  reduction  in  tlie  transverse  compressive  sti'oss 
necessary  to  produce  buckling. 

2.  An  interaction  formula  in  general  use  for  rectangular 
plates  in  combined  shear  and  compression  is  decidedly  con¬ 
servative  when  applied  to  an  infinitely  long  plate  in  slicar 
and  transverse  compression. 
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APPENDIX  A 

SYMBOLS 


Oi,  Ci  functions  of  edge  restraint  coefficient  «  given  in 
appendix  C 

D  fiexural  stiffness  of  plate  per  unit  length,  in-lb 


elastic  modulus  of  material,  psi 
compressive  force  per  unit  length,  Ib/in. 
shearing  force  per  unit  length  of  plate,  Ib/in. 
rotational  stiffness  per  inch  of  restraining  member 
at  edge  of  plate,  Ib/radian 

Tc  work  done  by  compressive  force  per  half  wave 

length,  in-lb 

T,  work  done  by  shear  force  per  half  wave  length, 

in-lb 

Vi  strain  energy  in  plate  per.  half  wave  length,  in-lb 

T7g  strain  energy  in  edge  restraint  per  half  wave  length, 

in-lb 

Y  function  of  y  associated  with  deflection  of  plate 

during  buckling 

b  width  of  plate,  in. 

bi  width  of  plate  in  oblique  coordinate  system  of  ref¬ 

erence  1,  in. 

c  function  of  a,  jS,  and  X 

ke,  k,  critical  compressive  and  shear-stress  coefficients, 
respectively 


m  root  of  a  characteristic  equation  of  appendix  B 

fitJitJi  functions  of  restraint  coefficient  e  given  in 
appendix  C 

t  thickness  of  plate,  in. 

vb 

w  displacement  of  buckled  plate  from  original  position 

Wo  amplitude  of  assumed  wave  form  of  buckle 

X  longitudinal  coordinate  of  plate 

y  transverse  coordinate  of  plate 

a,  jS  functions  of  X,  7,  and 

6  nondimensional  coefficient  of  edge  restraint 

7  one  of  two  parameters  determining  buckle  form 

X  one  of  two  parameters  determining  buckle  form 

(half  wave  length  of  buckle,  in.) 
ft  Poisson’s  ratio 

<r  direct  stress 

<r,  transverse  direct  stress,  psi 

T  shear  stress,  psi 

<j>  angle  between  buckle  node  and  y-axis 

0=tan  0 

Ri,  Rt,  Rt  stress  ratios 
R,  shear-stress  ratio 

Re  longitudinal  direct-stress  ratio 

er  critical  (used  as  subscript) 
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APPENDIX  B 

SOLUTION  BY  DIFFERENTIAL  EQUATION 


Statement  of  problem. — The  exact  solution  for  the  critical 
stress  at  which  buckling  occurs  in  a  flat  rectangular  plate 
subjected  to  combined  shear  and  compression  in  its  own 
plane  may  be  obtained  by  solving  the  differential  equation 
that  expresses  the  equilibrium  of  the  buckled  plate.  The 
plate  is  assumed  to  be  infinitely  long,  and  equal  elastic 
restraints  against  rotation  are  assumed  to  be  present  along 
the  two  edges  of  the  plate. 

Differential  equation. — ^Figure  6  shows  the  coordinate 
system  used.  The  differential  equation  for  equilibrium  of  a 
flat  plate  imder  shear  and  transverse  direct  stress  is  (from 
reference  7) 

It  is  convenient  to  write  o-„  and  t  in  terms  of  the  dimension¬ 
less  buckling  coefficients  k^  and  k,  by  means  of  the  relations 


keir^D  \ 

''  bH 

"  m 


(B2) 


Substitution  of  the  expressions  for  <r,  and  t  from  equations 
(B2)  in  equation  (Bl)  gives 

5*10  ,  „  d*to  ,  ,  2v*h, 

Solution  of  differential  equation. — If  the  plate  is  infinitely 
long  in  the  a:-direction,  aU  displacements  must  be  periodic 
in  X  and  the  deflection  surface  may  be  taken  in  the  form 


w=Ye  ^ 


(B4) 


where  Fis  a  fimction  of  y  only  and  X  is  the  half  wave  length 
of  the  buckles  in  the  aniirection. 
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Figcex  5. — C  oordinate  system  used  in  appendix  B  for  an  Infinitely  long  plate  under  combined 
shear  and  transverse  compression. 


Substitution  of  the  expression  for  w  from  equation  (B4) 
in  the  differential  equation  (B3)  gives  the  following  as  the 
equation  that  determines  F: 

d^Y  , /'x^kc  2v^(PY  ,  2ic^ik,dY  ^ 

A  solution  of  equation  (B6)  is 

F=e  ’> 

where  m  is  a  root  of  the  characteristic  equation 

m*-b[2  (y  J-  CB6) 

Except  for  the  substitution  of  2  2  ’ 

equation  (B6)  is  identical  with  equation  (A-6)  of  appendix 
A  of  reference  8,  in  which  equation  (B3)  of  this  appendix 
was  solved  with  ^*=0.  With  this  change,  aU  ,4he  results 
obtained  in  appendix  A  of  reference  8  are  applicable  here. 
The  stability  criterion  for  combined  compression  and  shear 
is  therefore  the  same  in  form  as  that  for  shear  alone,  given 
by  equation  (A-19)  of  reference  8,  which  is 

2aiS^47*— ^  (cosfr  2a  cos  2j8— cos  47)— j^47*(j8*— ce*) 

—  (|8*+a®)*— (47*— /3*-t-a*)  ^  sinh  2a  sin  2j8 

+ela(47*-|-a*-|-j3*)  cosh  2a  sin  2/S 

-t-iS(47*— a*— j8*)  sinh  2q:  cos  2^— 4aj37  sin  47]=0  (B7 

The  relation  between  k,  and  a,  /3,  and  7  is  also  the  same  in 
form  as  that  in  equation  (A-23)  of  reference  8,  namely. 


.  87(0*+^*) 


(B8) 


In  the  present  report,  however,  a  and  /3  have  the  following 
values: 


where 


,  ,  1  / xb\  X*  , 

®“'^+4(,xy  8** 


(B9) 
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As  in  reference  8,  the  restraint  coefficient  «  is  defined  herein 
by  the  relation 

Sb 

where  S  is  the  ratio  of  a  sinusoidally  applied  moment  to  the 
resulting  sinusoidally  distributed  rotation  of  the  restraining 
element  measured  in  radians. 

Evaluation  of  k,  corresponding  to  a  selected  value  of  kc. — 
The  procedure  for  evaluating  k„  after  values  of  ke  and  e 
have  been  chosen,  is  as  follows:  A  value  of  bj\  is  selected;  a 
series  of  values  of  y  are  assumed  until  one  is  found  that, 
together  with  the  corresponding  values  of  a  and  /3  computed 
from  equation  (B9) ,  satisfies  equation  (B7) ;  k,  is  then  com¬ 
puted  from  equation  (B8),  Another  value  of  bj\  is  selected; 
a  new  set  of  values  of  y,  a,  and  P  is  found  that  satisfies  the 
stability  criterion;  and  a  new  value  of  k,  is  computed.  The 
entire  process  is  repeated  until  the  minimum  value  of  k, 
can  be  found  from  a  plot  of  k,  against  6/X.  When  «  is  a 
function  of  6/X,  e  must  be  reevaluated  each  time  a  different^ 
value  of  6/X  is  selected.  The  minimum  value  of  k,  and  the 
chosen  value  of  ke,  when  inserted  in  equations  (B2),  give  a 
critical  combination  of  shear  and  direct  stress. 

Evaluation  of  k,  when  ke  has  value  corresponding  to 
buckling  as  Euler  strip. — One  critical  combination  of  shear 
and  compression  is  simply  A,=0  and  ke  equals  the  value 
corresponding  to  budding  as  an  Euler  strip.  The  curves 
giving  critical  stress  combinations,  however,  did  not  appear 
to  be  approaching  this  point  as  their  construction  pro¬ 
gressed.  It  was  therefore  necessary  to  determine  whether 
values  of  k,  other  than  zero  are  critical  when  the  Euler  com¬ 
pressive  stress  is  reached.  The  determination  of  k,  when 
ke  reaches  the  value  at  which  the  plate  buckles  as  an  Euler 
strip  requires  special  treatment,  because  k„  given  by  equa¬ 
tion  (B8),  becomes  indeterminate  when  the  wave  length  be¬ 
comes  infinite  as  su^estcd  by  the  energy  solution.  The 


result  that  X  becomes  infinite  when  ke  takes  its  Euler  value 
is  readily  checked  from  equation  (Cll)  for  the  special  case 

of  e=0;  for  this  case  p=2=r=^  and  (^:e)«r=l.  From 

equation  (B8)  it  is  clear  that,  if  is  to  remain  finite  when 
the  wave  length  approaches  infinity,  either  7->0,  case  (a),  or 
a*-!- case  (b). 

For  case  (a),  when  «=0,  it  follows  from  equations  (B9) 
that,  to  small  quantities  of  the  second  order. 


where 


(BIO) 


If  the  values  of  a  and  |3  from  equations  (BlO)  are  substi¬ 
tuted  in  equation  (B7)  and  the  resulting  equation  is  e.x- 
panded,  with  only  the  lowest  powers  of  u  and  y  retained, 


,  4x*u® 

64-  67r* 


(BID 


Substitution  of  the  values  of  a  and  jS  from  equations  (BlO) 
and  (Bil)  in  equation  (B8)  gives  as  the  final  result  for  «=0 

f.-+  _ 

*  *V64-6x* 

Case'IB)  can  be  analyzed  by  a  similar  method,  but  the 
analysis  is  quite  complicated  because  terms  of  third  order 
must  be  retained.  For  €=0  and  e=oo,  case  (a)  and  case 
(b)  were  found  to  lead  to  e.xactly  the  same  result  for  k,. 
A  value  of  k,  other  than  zero  when  ke  takes  its  Euler  value 
may  be  found  in  the  same  manner  for  other  values  of  edge 
restraint..  For  any  value  of  the  restraint  coefficient  e, 


ke* 

_(etVF.+2vVFc-H^^j-^ 

1  sinTV^+ 

1 

1 

2  " 
1  cos  (cos  xV^e— 1) 

1 

Isinx-V 

r.+( 

3x*  .ePke  3  D 
4  ‘  ‘  6  4  J 

|COS  xV^'c+l 

('4-b8E+2E*-t- 

^^(cosxVF.-l) 

(B12) 


where  ke  has  the  value  corresponding  to  Euler  strip  buckling, 
at  this  value  of  e.  The  relationship  between  this  Euler  value 
of  kc  and  the  corresponding  e  is  given  by  the  equation 
(from  reference  9) 


—  vJFe 


When  ke  reaches  the  Euler  value,  the  critical  shear-stress 
coefficient  k,  can  tlierefore  be  either  0  or  the  value  given  by 
equation  (B12).  The  conclusion  that  k,  can  also  have  any 
value  between  these  limits  is  plausible  on  the  basis  of  the 
following  physical  considerations.  The  shear  stress  does  no 


work  during  buckling  when  the  stress  condition  is  such  that 
the  plate  buckles  with  an  infinite  wave  length.  TIio  effect 
of  shear,  furthermore,  is  to  reduce  the  wave  length  to  a  value 
of  the  order  of  the  width  of  the  plate.  The  wave  length  at 
the  time  buckling  occurs  is  infinite,  however,  when  the  plate 
is  either  in  pure  compression  or  at  the  value  of  k,  satisfying 
equation  (Bl2).  This  fact  means  that,  for  values  of  k, 
between  0  and  that  given  by  equation  (Bl2),  the  shear 
stress  is  not  great  enough  to  force  buckling  in  short  waves 
and  therefore  does  not  assist  in  producing  buckles.  In  this 
rai^e  of  shear  stress  the  compressive  stress  necessary  to 
produce  bucklii^  is,  consequently,  the  Euler  stress. 


APPENDIX  C 


SOLUTIOX  BY  ENERGY  METHOD  FOR  EDGE  RESTRAINT  INDEPENDENT  OF  WAVE  LENGTH 


The  critical  stress  is  determined  on  the  basis  of  the  prin¬ 
ciple  that  the  elastic-strain  enei^  stored  in  a  structure 
during  buckling  is  equal  to  the  work  done  by  the  applied 
loads  during  buckling.  If  the  structure  under  consideration 
is  an  infinitely  long  plate  under  combined  shear  and  edge 
compression  with  edges  elastically  restrained  against  rota¬ 
tion,  this  equality  may  be  mitten 


Tc+T.=  V,+  V,  (Cl) 

In  reference  1  an  energy  solution  was  given  for  the  type 
of  loading  shown  in  figure  1  (b).  The  deflection  function 
used  in  reference  1  is  also  suitable  for  application  to  the 
solution  of  the  t3rpe  of  loading  shown  in  figure  1  (a),  which 
is  the  loading  considered  in  the  present  paper.  The  values 
for  T„  Fi,  and  Fj  may  accordingly  be  taken  directly  from 
reference  1,  but  Tc  must  be  recomputed  to  apply  to  the  case 
of  transverse  compressive  stress. 

The  following  substitutions  are  used  to  transform  the 
energy  expressions  from  the  oblique  coordinates  of  reference 
1,  the  left-hand  terms  of  the  equations,  to  the  rectangular 
coordmates  (fig.  6)  used  in  the  present  paper,  the  right-hand 
terms  of  tlie  equations: 

V  y\ 

^rh\  (C2) 

cos  4>—b) 

For  brevity  the  following  notation  is  also  adopted 

tan  (j>=d 
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Fiol'BE  6.— Coordinate  ssrstem  used  and  wave  form  assumed  for  energysolutkm  In  appendix  C ; 
Indlned  lines  Indicate  nodal  lines  of  buckles. 


By  use  of  equations  (C2),  the  expressions  from  reference  1 
that  are  used  in  the  present  paper  may  be  rewritten  as 
follows : 


T,=Wo^ 


Nxi/ic^bd  j. 
2\ 


(C3) 


yi=W  [0’(l  +  fi*)Vi+2(l  +  3fl*)/*+0J  /a]  (04) 


Fa='u;o* 


(C5) 


where 

and  e  is  the  restraint  coefficient  defined  in  appendix  B. 

The  work  done  by  the  compressive  force  per  half  wave 
length  may  be  written 

2  2 


As  in  reference  1,  the  assumed  deflection  function  is 

W=Wo  +f)  T J  f 

When  the  expre.^sion  for  w  from  equation  (C7)  is  substi¬ 
tuted  in  equation  (C6)  and  the  indicated  operations  are  per¬ 
formed, 

Tc=w^^(jJ,+6>^)  (C8) 

Values  from  equations  (C3),  (C4),  (C5),  and  (C8)  are  now 
substituted  in  the  buckling  equation  (Cl).  The  use  of  the 
equations 

N  —k 

iViy - K, 

N  —k  ^ 

XVy - Ac 


eliminates  Nxv  and  N^.  The  resulting  equation  gives  the 
critical  combination  of  stresses  and  may  be  written  as 

(i+TO+p/8+2  (1+30/* 

+Pp-*o(p/*+0‘i)]  (C9) 

This  equation  shows  that,  for  a  selected  value  of  the  criti¬ 
cal  shear  stress  depends  upon  the  wave  length  and  the  angle 
of  the  buckle.  Since  a  structure  buckles  at  the  low'est  stress 
at  which  instability  can  occur,  k,  is  minimized  with  respect 
to  wave  length  and  angle  of  buckle.  The  minimum  value  of 
k,  with  respect  to  value  of  wave  length  is  determined  from 
the  condition 

^■-0  (CIO) 
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which  gives  (when  e  does  not  depend  on  wave  length) 


(Cll) 


Substitution  of  this  value  of  wave  length  in  equation  (C9) 
gives 


(C12) 


The  minimum  value  of  k,  with  respect  to  angle  of  buckle  is 
found  from  the  condition 


^•=0 
be  ^ 


which  gives 


\—kJi 

l/s 

(  1 

A 

(C13) 


(C14) 


If  this  value  of  6  is  substituted  in  equation  (Cl2),  the  final 
result  is  the  following  interaction  formula,  in  which  k,  is 
given  in  terms  of  ke  and  the  edge  restraint  e: 


fe,*=4C?i+30,  {Ct-kc)-\-iAC,-kc)^!2(2C^=^^  (Cl5) 


where 


and 
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